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Abstract 



We investigate the possibility that the antiproton-to-proton flux ra- 
tio which is measured in cosmic rays may be generated by neutralino- 
neutralino annihilation in the galactic halo. We study the most general 
compositions for the relic neutralinos. We compare our results with the 
present experimental sensitivity and find that the theoretical predictions 
are at the level of the current experimental limits for some regions in the 
parameter space of the model. We expect that the future measurements 
of the p/p will provide very useful information, complementary to the 
ones obtainable with other experimental means. 



1 Introduction 



One of the best motivated particle candidates for cold dark matter is provided 
by supersymmetryLl. The lightest supersymmetric particle (LSP), being stable 
provided R-parity is conserved, is, in a large region of the supersymmetric 
parameter space, the lightest neutralino (\)- This state is defined as the lowest- 
mass superposition of the gaugino (13, W 3 ) and the higgsino {H\,H.2) fields : 

X = a 1 B + a 2 W 3 + a z H 1 +a i H 2 . (1) 



The possibility of detecting neutralinos in our galactic halo has been studied 
extensively in many different waysa. Direct detection of the nuclear jrecoil 
induced by neutralino-nucleus elastic scattering is, both theoretically 13 and 
experimentallyQa, under investigation. The indirect detection of high-energy 
neutrinos coming from the centre of the Earth or the Sunt), or the indirect 
detection of exotic components in the primary cosmic rays, such as 7, e + and 



p, have been studied in detail Ira. Here we examine the latter possibility, and 
analyse the parameter space of the supersymmetric model taking care of the 

Sew experimental limits on supersymmetric particles coming from CERN LEP 
. We also include some recent progresSpQn the treatment of the propagation 
of cosmic antiprotons p in the Galaxy EJ, considering the properties of the 
diffusion model instead of using the standard leaky box approach. 



2 Theoretical framework 



We recall in this section only the main features of the model to which we 
refer: the Minimal Supersymmetric Standard Model (MSSM). For a detailed 
description of the model and for a discussion on the range of variation of the 
various parameters that enters in the current analysis see: A. Bottino, these 
Proceedings. 

The neutralino mass and the coefficients at of Eq.(Q) depend on the pa- 
rameters: /i (Higgs mixing parameter), Mi, M2 (masses of B and of W 3 , 
respectively) and tan/3 = v u /vd (v u and Vd are the v.e.v.'s which give masses 
to up-type and down-type quarks). As for the masses of the particles enter- 
ing in the amplitudes for the neutralino-neutralino annihilation and for the 
neutralino-quark scattering (Higgs bosons, sfcrmions), we do not make use of 
theoretical assumptions implied by supergravity schemes, but we only take 
into account constraints due to present experimental limits. In the case of 
the neutral Higgs bosons (the two CP-even bosons: h, H (of masses rrih, run 
with run > nrih) and the CP-odd one: A (of mass vtia)) we take-m^ as a 
free parameter only bounded by the present experimental results Ej. As for 
the sfermion masses and in order to deal with a relatively small number of 
parameters, we consider a common soft scalar mass mo as a free parameter 
(with the possible exception of the top squark). Another quantity which we 
take as a free parameter is the top trilincar coupling A t . Moreover we assume 
the standard GUT relation among the gaugino masses 

Mi : Mi : M3 — ot\ : ct2 ■ «3 (2) 

which implies at raz scale that M\ ~ O.5M2. 

Once defined our parameter space, we impose the experimental constraints 
deduced from the new particles searches at accelerators, including the new 
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lower limits on the chargino mass obtained at LEPIlEI, and the measurement 
of the rare decay b — > s-yliS from CLEO. 

hen compute the neutralino relic abundance in the standard way 
and we exclude from the parameter space those configurations that 
do not fulfil the cosmological bounds ilh 2 < 1. 




2.1 Relic density 



The neutralino relic abundance fl x h 2 is evaluated following the standard pro- 
cedure, which gives essentially Vt x h 2 oc< cr ann w >j~ t , where < a ann v >; nt is the 
thermally-averaged annihilation cross section, integrated from the freeze-out 
temperature to the present temperature. The standard expansion < a^nV >= 
a + bx + ... may be employed, with x — T/m x , except at s-channel resonances 
(Z, A, H, h) or at the opening of a new final state chann£l, where a more pre- 
cise treatment has to be used for the thermal averages. In the evaluation 
of < GannV > the full set of annihilation final states (// pairs, gauge-boson 
pairs, Higgs-boson pairs and Higgs-gauge boson pairs), as well as the complete 
set of Born diagrams, are taken into account Ej. We note that the constraint 
^x^ 2 — 1 can be very effective in cutting the parameter space especially for 
small and intermediate values of tan (3, but is not really restrictive for large 
values of tan (3. 



3 The antiproton spectrum 



Let us turn now to the evaluation of the antiproton signal due to neutralino 
annihilation in the galactic halo. The differential rate (per unit volume and 
per second) for the production of antiprotons from x~X annihilations is given 
by 

S(E P ) = (av)f(E p ) (£-) , (3) 

in which Ep denotes the antiproton energy, a is the x~X annihilation cross 
section and v is the neutralino velocity in the galactic halo. The neutralino 
density p x is a function of the position r in the galactic halo. For a single 
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annihilation, the antiproton energy spectrum is 

1 da( XX ^p + X) v {F) ( dflA 

f ^ = a dEp -jrj xf [d*) ' ( ) 

where F describes the x~X annihilation final state and B ^ is the branching 
ratio into the quarks or gluons / in the channel F. The differential distribution 
of the antiprotons generated by the hadronization of quarks (with the exception 
of the top quark) and of gluons, is denoted by dNp /dEp and depends on the 
nature of the species /. In Eq.(jj|), (av) and f{Ep) depend on the neutralino 
properties. The antiproton production rate also depends on the distribution 
p x of neutralinos inside the galactic halo. 

For the p differential distribution f(Ep), we have evaluated the branching 

ratios B ^ for all annihilation final states that may produce antiprotons, i.e. 
direct production of quarks and gluons, generation of quarks through the in- 
termediate production of Higgs bosons, gauge bosons and the top quark. The 
distributions dNp /dEp from the hadronization of quarks (with the exception 
of the top quark) .and gluons have been computed by using the Monte Carlo 
code JETSET 7.2B. 

The neutralino halo distribution p x is taken to be spherically symmetric and 
is given by the standard expression: 

Px(T) = Px(®) ^r^f , (5) 

where a = 3 kpc is the core radius of the dark matter halo. Particular care 
must be taken about the local neutralino density p x (&), which depends on 
on the distribution of dark matter in the galactic halo as well as on the LSP 
properties. On the one hand the presence of a clumpy structure over the 
smooth distribution || may well enhance the expected signalEix On the other 
hand if the big-bang relic density fl x h 2 , which we evaluate following the method 
previously discussed, is too small to account for the cosmological dark matter, 
the density of neutralinos in the galactic halo should be corrected by a factor of 
£. The latter deals with the fact that the neutralino density is less than the halo 
density whenever fl x h 2 is smaller than a minimal value of, say, (flh 2 ) m i n — 
0.03, which is compatible with the observed rotation curve of the Galaxy. At 
cosmological distances, this ratio is given by 

^x h2 

s ~ {<nh 2 ) 

nini 
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Figure 1: The ratio p/p as a function of the neutralino mass for the simple scan of the param- 
eter space asrdpscribed in the text. The horizontal line represent the current experimental 
upper bound E9. 

The antiproton flux expected at Earth can be written as 



*p(Ep) = ^<<rv>f(Ep)(£^ v P r p (E p ) (7) 

where v p and Tp(Ep) are the p velocity and confinement time. 

As far as Tp(Ep) is concerned we employ the results of Ref0 and we in- 
troduce an energy dependence of Tp as deduced from a detailed calculation a 
diffusion model to treat the p propagation. 

To compare predictions with experimental data we have to take into ac- 
count the solar modulation of the antiprotons as welLas the modulation of the 
primary protons. We have used the results of Ref. Ej changing the values of 
the parameters of. the modulation model in order to properly compare with 
the experiments E9. 
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Figure 2: The ratio p/p as a function of the neutralino relic density Q, x h? for the simple scan 
of the parameter space as piascribed in the text. The horizontal line represent the current 
experimental upper boundc3. 



4 Results 



In Fig.l we display the ratio p/p as a function of the neutralino mass for a 
simple scan of the parameter space, where tan/3 = 1.1 and tan/3 = 55, and 
itia — 500 GeV and tua — 50 — 60 GeV. This choice represents almost the two 
boundaries for the range of variation of these parameters. The horizontal line 
represents the present level of sensitivity in the-energy range 0.25 GeV < T < 1 
GeV as measured by the IMAX collaboration^], that is 

{^j < 7.1 x 10" 5 95 % C.L. . (8) 

We notice that for some configurations of the model we are at the level of the 
experimental result. 

In Fig. 2 we show the p/p flux versus the cosmological relic density fl x h 2 . 
This plot illustrates the dependence of the signal on the rescaling procedure 
we have adopted, here we recall that the value (Qh 2 ) min — 0.03 has been used. 
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Figure 3: The ratio p/p as a function of the kinetic energy T. The solid line represent four 
different values of the % mass (from left to right: 40 GeV, 50 GeV, 60 GeV, 120 p^V). The 
dashed line represent the expected background due to secondary production of plH 



In Fig. 3 we display the p/p ratio as a function of the kinetic energy T. The 
different curves represent four different neutralino masses and compositions. 
In the same plot is shown, as a dotted line, the expected p/p background of 
the secondary production by spallation of the primary cosmic ray flux (taken 
fromRef. E3). 



A possible clear signature to discriminate signal over background is the 
energy dependence of the two fluxes on T: the former being a milder function 
of T at low energies with respect to the latter, as it is clearly seen in Fig. 3. 
This methods will be useful at the forthcoming generation of experiments such 
as AMS EJ or PAMELA E3 that will collect a much larger number of p with 
respect to the present experimental situation. 



In summary we have shown that the production of p's from \~X annihila- 
tion in our galactic halo is, in some regions of the MSSM parameter space, at 
the level of the present experimental sensitivity. The forthcoming generation 
of detectors will be able in the near future to investigate large regions of the 
MSSM by discriminating the signal over the background with an analysis of 



7 



the energy dependence of the measured p/p. 
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